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Background. Gene expression profiling test scores have primarily been used to identify heart transplant recipients 
who have a low probability of rejection at the time of surveillance testing. We hypothesized that the variability of gene 
expression profiling test scores within a patient may predict risk of future events of allograft dysfunction or death. 
Method. Patients from the IMAGE study with rejection surveillance gene expression profiling tests performed at 1- to 
6-month intervals were selected for this cohort study. Gene expression profiling score variability was defined as the 
standard deviation of an individual's cumulative test scores. Gene expression profiling ordinal score (range, 0-39), 
threshold score (binary value=l if ordinal score >34), and score variability were studied in multivariate Cox regression 
models to predict future clinical events. 

Results. Race, age at time of transplantation, and time posttransplantation were significantly associated with future 
events in the univariate analysis. In the multivariate analyses, gene expression profiling score variability, but not 
ordinal scores or scores over threshold, was independently associated with future clinical events. The regression co- 
efficient P values were <0.001, 0.46, and 0.773, for gene expression profiling variability, ordinal, and threshold scores, 
respectively. The hazard ratio for a 1 unit increase in variability was 1.76 (95% CI, 1.4-2.3). 

Discussion. The variability of a heart recipient's gene expression profiling test scores over time may provide prog- 
nostic utility. This information is independent of the probability of acute cellular rejection at the time of testing that is 
rendered from a single ordinal gene-expression profiling test score. 
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The standard care of adult heart transplant recipients has 
included surveillance for rejection with serial endo- 
myocardial biopsies (1-3). Because of the risks and discom- 
forts associated with the biopsy procedure, a noninvasive gene 
expression profiling test (AlloMap) of peripheral blood was 
developed to identify heart transplant recipients who have a 
low probabihty of rejection at the time of surveillance testing 
(4-6). We have noticed an association of low variability of 
gene expression profiling scores within an individual and the 
clinical stability of the individual ( 7) . Therefore, we hypoth- 
esized that the variability of gene expression profiling scores 
within individuals may predict risk of future allograft events. 

To test our hypothesis, we utilized the IMAGE (Inva- 
sive Monitoring Attenuation by Gene Expression Profiling) 
study database, which provided independently adjudicated 
clinical outcome events observed during surveillance of 602 
heart transplant recipients (8). In IMAGE, recipients at least 
6 months posttransplantation were randomized 1:1 to either 
surveillance with routine biopsy or gene expression profiling 
testing. The patients assigned to gene expression profiling 
underwent a surveillance biopsy only if the gene expression 
profiling score was equal to or above the specified threshold 
of 34. A primary outcome event was defined as an episode of 
rejection with hemodynamic compromise, graft dysfunction 
due to other causes, death, or retransplantation. Over a median 



of 19 months of observation, 297 patients monitored with gene 
expression profiling and 305 patients monitored with biop- 
sies had similar 2-year cumulative rates of events (14.5% 
and 15.3%, respectively; hazard ratio with gene expression 
pro filing= 1.04; 95% confidence interval, 0.67-1.68) (8). 

RESULTS 

The demographics of the study cohort (N=369) were 
representative of the features of the parent IMAGE population 
(n=602) (Table 1). The cohort was predominandy white (80%) 
men (83%), 54 years of age at the time of transplantation. Study 
related surveillance began after 12 months posttransplantation 
in 86% of the cohort. 

Within Patient Gene-Expression Profile Scores 
and Variability 

The cohort had a mean of 4.4 (SD, +1.7) gene expression 
profiling tests and a mean ordinal score of 30.9(SD, ±0.7) (see 
Table SI and Figure SI, SDC, http://links.lww.com/TP/A924). 
The gene expression profiling variability was 1.0 (SD, +0.8); 74 
patients (20%) had variability less than 0.5, 160 patients 
(43%) had variability greater than 1.0, 65 patients (18%) had 
variability greater than 1.5, and 26 patients (7%) had vari- 
ability greater than 2.0. The variability was 1.6 (SD,+1.4) units 
for the event and 1.0 (SD, +0.7) units for the no events 



TABLE 1 . Demographics of study cohort compared to whole IMAGE population 






IMAGE population 


Multivariate analysis 


Characteristics 


(N=602) 


study cohort* (N=369) 


Age, yr: mean+SD 


54.1+12.8 


54+12.9 


Male sex, no. (%) 


493 (82%) 


305 (83%) 


Race or ethnic group, no. White (%) 


468 (78%) 


294 (80%) 


Indication for cardiac transplantation, no. (%) 






Coronary artery disease 


257 (43%) 


160 (43%) 


Nonischemic cardiomyopathy 


307 (51%) 


185 (50%) 


Other 


38 (6%) 


24 (7%) 


Interval between transplantation and randomization, no. (%) 






6-12 mo 


88 (15%) 


52 (14%) 


13-36 mo 


412 (68%) 


264 (72%) 


37-60 mo 


101 (17%) 


52 (14%) 


Cytomegalovirus status, no. (%) 






Donor and recipient positive 


237 (39%) 


138 (37%) 


Donor and recipient negative 


91 (15%) 


47 (13%) 


Donor positive and recipient negative 


137 (23%) 


98 (27%) 


Donor negative and recipient positive 


108 (18%) 


77 (21%) 


Unknown 


29 (5%) 


9 (2%) 


Use of ventricular assist device before transplantation, no. (%) 


115 (19%) 


69 (19%) 


Induction therapy (any), no. (%) 


352 (58%) 


220 (60%) 


Immunosuppressive therapy, no. (%) 






Cyclosporine 


162 (27%) 


94 (25%) 


Tacrolimus 


436 (72%) 


280 (76%) 


Mycophenolate mofetil or mycophenolic acid 


486 (81%) 


306 (83%) 


Prednisone 


271 (45%) 


168 (46%) 


Left ventricular ejection fraction at first study visit 






MeaniSD 


60.4+6.3 


60.6+6.4 


No. of events, no. (%) 


67 (11%) 


30 (8%) 
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TABLE 2. Univariate Cox proportional hazards regression analyses of features associated with clinical events" 
Independent variable Regression coefficient P Hazard ratio (95% CI)'' 



Sex 

Race: nonwhite 

Age at time of transplant (log year) 



Time posttransplant (log days) 



No. rejection episodes in the first year 
»Cytomegalovirus (CMV) serologic status 

(ordinal: 0 to 3: D-/R-, D-/R+, D+/R+,D+/R-.)^ 



-0.07 
1.05 
-2.57 



-2.39 



0.049 
0.09 



0.83 
<0.001 
<0.001 



0.005 



0.72 
0.48 



0.9 (0.5-1.7) 
2.9 (1.8-4.7) 
0.64 (0.50-0.82) for a 30-year-old 

versus a 20-year -old 
0.84 (0.77-0.93) for a 70-year-old 
versus age 60) 

0.69 (0.53-0.89) 
for (12 mo vs. 18 mo) 
0.85 (0.75-0.95) for (30 mo vs. 36 mo) 
1.1 (0.8-1.4) 
1.1 (0.9-1.4) 



" Clinical events of rejection with hemodynamic compromise, graft dysfunction due to other causes, death, or retransplantation. 
Hazard ratio: values less than 1 indicate risk decreases as variable increases; the associated regression coefficient is negative. 

0=D-/R- is donor and recipient negative, 1=D-/R+ is donor negative and recipient positive, 2=D+/R+ is donor and recipient positive,3=D+/R- is donor 
positive and recipient negative. 

Note: Each row in Table 2 corresponds to a separate Cox proportional hazards regression analysis. 



subgroups, respectively (P=0.0126, Mann- Whitney test). The 
within-patient variability scores were weakly associated with 
the respective within-patient median ordinal scores (Spearman's 
correlation coefficient, r=0.14, _P<0.001). 



Base Clinical Risk Factors Model 

The univariate analyses revealed that nonwhite race, 
younger age at time of transplantation, and earlier time 
posttransplantation at study entry were significantly asso- 
ciated with a higher risk of events (Table 2). A nonwhite was 
2.9 times more likely to have an event than a white. Sex, 
number of rejection episodes experienced in the first year 
posttransplantation, and cytomegalovirus serologic status 
were not significantly associated with events. The significant 
risk factors were included in the multivariate models below. 



Multivariate Models of Gene Expression 
Profiling Test Information to Predict Future 
Clinical Events 

Among the multivariate models, only gene expression 
score variability was significantly predictive of future clinical 
events (regression coefficient, P<0.001; Table 3). Gene ex- 
pression score variability had the lowest estimated informa- 
tion loss (Akaike information criterion (9)) and a predictive 
accuracy (concordance index) of 0.69 (95% CI, 0.61-0.76). 
Gene expression score variability was significant even when 
controlling for gene expression ordinal score (_P<0.001; see 
Table 82, SDC, http://links.lww.com/TP/A924. The C-index 
for the base clinical risk factors was 0.68 (95% CI, 0.59-0.78). 

The potential clinical utility of the variability score is il- 
lustrated in four cases fi-om the IMAGE study (Fig. 1). The 
hazard ratio associated with the variability score and adjusted 



TABLE 3. Multivariate models to predict future clinical events'" 

Cox proportional hazards regression results of surveillance information added to base clinical model 




Regression 




Hazard ratio 




C-Index 




coefficient (SE) 


P 


(95% CI)*^ 


AIC 


(95% CI) 


Base clinical model'' 








329.5 


0.68 (0.59-0.78) 


Gene expression profiling test score variability: 


0.57 (0.12) 


<0.001 


1.76 (1.38-2.25) 


318.0 


0.69 (0.61-0.76) 


(standard deviation of all previous 












gene expression profiling ordinal 












scores since study entry)'' 












Gene expression profiling test 


-0.078 (0.10) 


0.43 


0.92 (0.76-1.12) 


330.9 


0.68 (0.59-0.77) 


ordinal score (0-39) 












Gene expression profiling 


0.12 (0.43) 


0.77 


1.13 (0.48-2.65) 


331.5 


0.68 (0.59-0.78) 



test threshold score >34 



" Clinical events of rejection with hemodynamic compromise, graft dysfunction due to other causes, death, or retransplantation. 

'' Race, younger age at time of transplantation, and earlier time posttransplantation were the variables included in the Base Clinical Model. 

' Hazard ratio: values less than 1 indicate risk decreases as variable increases; the associated regression coefficient is negative. 

Uses transformed gene expression profiling scores. If x is the AlloMap gene expression profiling ordinal score, we use y=2.451(log(x/(40-x))— 0.234). 
Note: Each row in Table 3 corresponds to a separate Cox proportional hazards regression analysis. 
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FIGURE 1. Examples of using gene expression profiling test score variability to predict future clinical events*. In each 
panel, the top inset graph shows the longitudinal visits and associated gene expression profiling ordinal test scores (the time 
of the first surveillance gene expression profiling test is treated as day 0). The index visit number and the associated 
variability score is indicated with the arrow. In the lower inset in each panel, the "probability of event" is plotted, along with 
95% CI (dashed lines). The risk (e.g., 2.6%) of an event in the next 12 months for the patient is shown with the red lines. 



for base clinical risk factors was applied to estimate the prob- 
ability of future events. The variability score is time varying: for 
example, in patient B, the variability score is 1.17 at visit 3 and 
1.10 at visit 4. Only one example variability score is applied to 
project the probability of an event in each illustration. 



Patient A 

This white man, transplanted at age 57, was 15 months 
posttransplant at the time of study entry. Patient A had 
relatively high ordinal scores (35, 35, 36, 36, and 35); be- 
cause the scores exceeded the nominal threshold (>34), each 
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triggered a per protocol endomyocardial biopsy. Patient A 
did not have biopsy-proven rejection in association with his 
34 or higher scores (biopsy grades of 0, 0, IR, 0 and IR, 
respectively (10)). At visit 5, his score variability (0.49) was 
at the 20th percentile of the study cohort. Based on the 
multivariate model, patient A with this low variability of 
gene expression profiling scores had a relatively low proba- 
bility (2.6 %) of a clinical event in the future 12 months 
from visit 5, although each of his individual ordinal scores 
exceeded the nominal threshold of 34, which indicated in- 
creased risk of moderate or severe (>grade 2R) acute cellular 
rejection. This case illustrates that a relatively higher prob- 
ability of acute cellular rejection associated with an ordinal 
score 34 or higher (i.e., a 6% to 10% positive predictive 
value, established from the CARGO population study) may 
coexist with a relatively low probability of future clinical 
events. The probability of future clinical events is adjusted 
for the patient's underlying clinical risk factors in this 
model. The risk for future events (defined clinically, without 
requirement of histologic evidence of rejection) based on 
the variability of an individual's gene expression scores may 
be independent and uncoupled from the risk of acute cel- 
lular rejection (defined by histology) based on single gene 
expression profiling ordinal scores at the time of testing. 

Patient B 

This white man, transplanted at age 34, was 15 months 
posttransplant at study entry. His ordinal scores (26, 23, 31, 
23, and 33) render a higher variability score (1.10 at visit 4, 
which is at the 64th percentile of the population) compared 
with patient A. Patient B's overall higher risk of having a 
future event (11.5%) after visit 4 is due to his younger age at 
time of transplantation (age 34) and his higher gene ex- 
pression profiling score variability compared with patient A. 
Patient B went on to incur a clinical event of hemodynamic 
compromise without rejection about 4 months after his visit 
4. Six months later ( — day 1050 posttransplant), he experi- 
enced another episode of hemodynamic compromise with 
histologic evidence of grade 3A acute cellular rejection. 

Patient C 

This African American woman, transplanted at age 57, 
was 16 months posttransplant at study entry. Her ordinal 
scores were 24, 26, 26, and 28. Her low score variability 
(0.46) at visit 4 was similar to that of patient A. The mul- 
tivariate model yields a net 12 -month probability of a future 
event of 7.8% forward from visit 4 for patient C. Patient C's 
race factor drives her higher risk for future events compared 
with patient A, the white middle-aged man with a similarly 
low score variability (0.49). In contrast, patient C has lower 
risk than patient B, a younger white man, with higher gene 
expression score variability. 

Patient D 

Patient D, a white man, was age 66 at the time of 
transplantation and 15 months posttransplant at study en- 
try. His ordinal scores were 29, 25, 31, 35, and 32. Based on 
only conventional clinical risk factors. Patient D would be 
considered to be at lower risk for future events compared 
with patient C, the African American woman. His one or- 
dinal score of 35 had an associated biopsy rejection grade of 
IR. However, from the multivariate model, the probability 



of a future event (7.1%) in patient D is similar to that of 
patient C (7.8%) and higher than patient A (2.6%). His gene 
expression profiling variability score was 1.34, which is at 
the 77th percentile of the study cohort. Patient D went on to 
incur a clinical event of rejection with hemodynamic com- 
promise, 3 months after his visit 5. 

DISCUSSION 

The above anecdotal cases help to illustrate the potential 
use of gene expression profiling variability to predict the 
probability of future clinical events in heart transplant re- 
cipients. For example, a recipient predicted to be at low risk 
for future events may become a candidate for further mini- 
mization of his/her immunosuppressive maintenance regi- 
men. Conversely, an individual predicted to be at higher risk 
for future events may receive further evaluation to detect 
possible underlying causes of the variability such as over- 
looked infections or noncompliance to medications. Our 
analyses so far have not enabled us to define the specific causes 
of a high variability of gene expression profiling scores. 

Future mechanistic studies, using a larger number of 
cases, should focus on identifying individual gene pathways 
associated with global score variability, to further understand 
the relationship between score variability and clinical out- 
comes. The expression levels of the 1 1 genes that are the core 
of the AJloMap gene expression profiUng test were chosen for 
their performance in linear discriminant analysis for the 
probability of the presence or absence of acute cellular rejec- 
tion (4). The changing expression levels of the genes over time 
in certain circumstances may reflect a change in the status of 
the recipient's immune system that is not directly associa- 
ted with acute cellular rejection. Up until now, the best cha- 
racterized specific factors that influence background gene 
expression profiling scores independent of acute cellular re- 
jection include corticosteroid doses above 20 mg/d of pred- 
nisone, which is known to suppress the ordinal score, and 
cytomegalovirus seropositivity of donor and/or recipient, 
which raises the average ordinal score relative to cytomega- 
lovirus seronegative recipients (11, 12). Race, sex, age, and 
blood levels of calcineurin inhibitors have relatively minor 
influence on AUoMap scores (13). 

In the fiiture, we may be able to identify further corre- 
lations of the gene expression profiling score, its subcomponent 
genes, and/or score variability to predict more specific events 
(e.g., morbidity or death fi^om cancer or cardiac allograft 
vasculopathy). It wiU be useful to further characterize gene ex- 
pression profiling score patterns in larger heart transplant co- 
horts that include more specific categories of clinical end point 
events collected over longer foUow-up periods to confirm and 
extend the current model findings. The next studies should also 
determine how variability is affected by the number and time 
intervals for repeat testing. The Cardiac Allograft Rejection 
Gene Expression Observational (CARGO II) observational trial 
and the early IMAGE (EIMAGE) wiU provide further data on 
the gene expression profiling score variability associated with 
earlier (2-6 months posttransplantation) and more frequent 
(interval of every 2—4 weeks) use of the gene expression 
profiling test; ClinicalTrials.gov Identifier NCT00761787 
and ClinicalTrials.gov Identifier NCT00962377, respectively 
(14-16). 
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The further investigation of gene expression profiling 
scores and score variabiHty on future cHnical events is a 
worthy endeavor because advances in the understanding of 
the determinants of long-term outcomes have been identi- 
fied to be among the highest priorities for future heart 
transplantation research (17). 

The small increase in C-index (0.01) when the vari- 
ability score is added to the base clinical model indicates 
limited improvement in discrimination of patients who will 
have events from those who will not have events. However, 
the 95 % confidential intervals range from 0.59 to 0.78, so it 
is possible that with a larger dataset, a larger difference in the 
C-index may be demonstrated. The variability factor seems 
to have promise, based on its highly significant regression 
coefficient and its favorable AIC score. The clinical utility of 
the variability score may be assessed based not only on its 
discrimination performance but also on its calibration per- 
formance (agreement between predicted and observed risks 
across subgroups with varying baseline risk) and/or the new 
test's reclassification performance (assessment of the pro- 
portion of individuals moved between risk categories rele- 
vant to treatment decisions) (18). These additional metrics 
are beyond the scope of this report. 

Our analysis does not distinguish if the variability 
observed is attributed to a temporal increasing (or decreas- 
ing) trend or to variability without a trend. The association 
of score trends with clinical outcomes has not been anec- 
dotaUy reported and was not proposed in the main clinical 
hypothesis of this paper. Trend analyses thus awaits assess- 
ment in a future study, where availability of more than four 
or five test results per patient would enhance the power of 
the analysis. 

Based on the results from the current study, we spec- 
ulate that independent from the established association of 
ordinal scores over a threshold value with higher risk of acute 
cellular rejection at the time of testing, the low variability of a 



set of gene-expression profiling scores within an individual 
may be associated with future clinical quiescence in an indi- 
vidual. A low gene expression profiling variability score may 
predict the clinical stability of an allograft recipient, inde- 
pendent of the fundamental interpretation of single ordinal 
test scores for risk of acute cellular rejection. 

In conclusion, the variability of gene expression profil- 
ing scores from an individual may help predict the risk of 
clinically defined future allograft dysfunction or death in the 
individual. This prognostic information is independent from 
and complementary to the original use of single gene ex- 
pression profiling ordinal test scores to estimate the proba- 
bility of histologically defined acute cellular rejection at the 
time of testing. To our knowledge, the current report is the 
first study of the utility of the variability of a biomarker score 
that is derived from a multivariable gene expression profile. 
Further development and validation of this prognostic test 
may enable improvement in the long-term management and 
clinical outcomes of heart transplant recipients. 

MATERIALS AND METHODS 
Study Design 

The current study was a cohort analysis of patients from the IMAGE study 
(8). SurveiUance visits occurred at 1- and 2-month intervals for patients who 
were between month 6 and 12 posttransplantation, and at 3-, 4- or 6-month 
intervals after the first year posttransplantation. The gene expression profil- 
ing test (AUoMap) was provided by XDx, Inc., Brisbane, CA (19). 

Patient Population and Surveillance Testing 

The current study cohort, n=369, was drawn from the 602 patients from 
13 US centers in the IMAGE study (8). These 369 patients were selected for 
inclusion because they had at least two gene expression profiling tests before 
an event or study end, and at least one endomyocardial biopsy and one 
echocardiogram were included (Fig. 2). During the course of surveillance of 
this cohort, primary clinical events were documented in 30 patients. Pa- 
tients in both IMAGE study arms provided blood samples for gene 



IMAGE Intent-to-Treat Cohort 
N = 602 



No Biopsy, No Echocardiogram LVEF, 
and/or <1 Gene-Expression Profiling 
Scores Prior to Event 
N = 15 



> 1 of Additional Surveillance Method 
Measurements While on Study 
N = 587 



> 1 Biopsy 
N = 433 



1 Ectiocardiograptiy LVEF 
N = 573 



Gene-Expression Profiling Scores 
N = 500 



> 1 LVEF, > 1 Biopsy, and >1 Gene-Expression Profiling Scores 
Included in Muitivariate Analysis 
N = 369 



FIGURE 2. Source of N=369 patients included in the modeling of prediction of future events of allograft failure. 
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expression profiling testing at surveillance visits; samples fi-om both arms of 
the IMAGE study were included in the current analyses. 

Clinical Events 

The definition of clinical events was taken directly from the IMAGE 
study. The composite end point was composed of a patient's first event(s) of 
the following: rejection with hemodynamic compromise, graft dysfunction 
due to other causes, death, and/or retransplantation. An independent end- 
points committee had adjudicated aU primary events (8). 

Statistics 

The association of clinical risk factors with clinical events was assessed 
in a univariate Cox regression model (R-2.13.1, R foundation, Vienna, 
Austria). The factors included patient's race, sex, age at time of transplan- 
tation, number of prior rejections in first year, cytomegalovirus serologic 
status, and time posttransplantation. Age at time of transplant and time 
posttransplantation were log transformed to better fit the risk of events in 
these analyses. These factors were selected based on prior published evi- 
dence of association with risk of poor allograft outcomes and availability of 
IMAGE data (8, 20). 

The factors that reached significance (P<0.05) in the univariate analyses 
were included with the gene expression profiling scores in the multivariate 
Cox proportional hazards models. These models provide a prediction of the 
risk of future events expressed as a hazard ratio from the gene expression 
profiling scores, adjusted for the base clinical risk factors. The underlying 
time variable in the Cox regression model was the time since enrollment in 
the study, rather than time since transplant. This analysis method was se- 
lected because, per the IMAGE study protocol inclusion criteria, no patients 
were enrolled until after 6 months posttransplantation. 

The gene expression profiling scores were analyzed in three ways: 1 ) the 
ordinal score, an integer ranging from 0 to 39; 2) the threshold score, 
assigned a binary value of 0 if the ordinal score was less than 34 and 
assigned a value of 1 if the ordinal score is 34 or higher; and 3) the vari- 
ability score, which is the standard deviation of aU ordinal scores measured 
within a patient, before incurring a clinical event. The method for com- 
puting the variability score is as follows. First, collect the gene expression 
profiling ordinal scores from the previous visits of an individual: for ex- 
ample, X =24, 23, 31, and 23. Second, transform each ordinal score by the 
inverse logit transformation, where for ordinal score x, y=2.451{log(x/(40-x))- 
0.234) is the transformed score. In the example, y=0.42, 0.17, 2.5, and 0.17. This 
transformation is necessary to better satisfy the modeling assumptions of the 
Cox model. Third, calculate variability score, which is the standard deviation 
of the transformed scores: 
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cr = . / y (x, — u) , where X: are the N=4 transformed scores, 
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and is the mean of the scores. In the example, the variability score is 1.11. 
The variability score is a time-varying covariate. Beginning with the second 
visit, one may compute the variability score and use that to predict the risk 
of events from then until the time of visit 3 or first event, whichever occurs 
first. If there is no prior event, then at visit 3, one may compute a new 
variability score based on three ordinal scores to predict an event up to the 
time of visit 4 (or first event) and so on. 

The multivariate models render a regression coefficient and hazard ratio for 
risk of future events for each of the gene expression profiling score variables, 
adjusted for the base clinical risk factors. The comparative performance of the 
three models was assessed by the Akaike information criterion (AIC) and the 
concordance index (C-index). The AIC is a measure of the estimated infor- 
mation loss in the statistical model; the smaller the AIC, the less is the infor- 
mation loss by the model relative to other models (9). The C-index is a 
generalization of the area under the receiver operator characteristic curve and 
measures how well the model performs in predicting fiiture clinical events (21). 
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